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ABSTRACT 


A comparison of two different autopilot configurations and their effect on missile 
response is presented. The comparison includes the steps taken in determining missile 
parameters from wind tunnel data and flight condition data. The missile parameters are 
coupled with two different autopilot configurations to determine any significant advan- 
tage of one configuration over the other. Pole placement is used in determining required 
autopilot feedback and feed forward gains. Simulations of each autopilot are conducted 


and the responses are conipared. 
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I. INTRODUCTION 


The purpose of an autopilot in a missile is to cause the missile to maintain stability 
as it travels along its flight path. Commonly the autopilot provides stability about the 
mussile’s roll, yaw and pitch axes. Autopilots generally use displacement gyros with 
feedback being provided by electronic sensors which detect axial displacernent due to 
nussile rotation about one or more of the above mentioned axes. 

Autopilots. as with most other man-made devices, can be constructed in a variety 
of different configurations. Some configurations may have advantages over others in 
terms of paraincters such as cost, weight or effectiveness, to name a few. 

In this research paper two pitch axis autopilot configurations are compared. The 
first 1s a three loop autopilot configuration, employing body rate feedback, lateral ac- 
celeration fecdback and an additional synthetic stability feedback loop as depicted in 


migure |. 


SYNTHETIC 
STABILITY 


LOOP 





Figure 1. Autopilot Configuration | 


The second is a two-loop autopilot configuration, employing proportional plus integral 


compensation in the feed forward path with body rate and lateral acceleration feedback 


as shown i [termes 





Figure 2. Autopilot Configuration 2 


In these two figures 7¢ 1s the commanded lateral acceleration in the pitch plane, 174 
is the measured lateral acceleration in the pitch plane, 6 is the missile control surface 
deflection angle and @ is the missile pitch angle rate. All the subscripted K variables are 
gains Which are constant at each design point. 

The purpose of the comparison of the two autopilot configurations is to determine 
what advantage, if any, one may have over the other in terms of missile response. 
Comparisons are made of the two autopilots for three different sets of flight conditions. 
Lach set of flight conditions is referred to as a design point. For example, design point 


one will refer to the first set of flight conditions. These flight conditions, or design 


bo 


points, are specified in Chapter Two. The desired autopilot response for each config- 
uration Is given as approximately 10 radians’second with a damping coefficient of 0.5. 

Pertinent missile parameters are derived in Chapter Two. using given wind tunnel 
data and flight condition data. Any needed parameters which cannot be derived or 
otherwise obtained through availible information, are assumed using sound engineering 
judgment. 

In Chapter Three the missile and autopilot transfer functions are developed. The 
feedback gains necessary to meet the desired autopilot design specifications are obtained 
allowing development of simulation studies. 

Simulations are developed and conducted using Dynamic Simulation Language 


(DSL) in Chapter Four and the results are analyzed. 


HW. DEVELOPMENT OF MISSILE PARAMETERS 


A. INTRODUCTION 

In order to compare the effects that the two autopilot configurations have on missile 
performance, determination of certain missile parameters must first be made. Specil- 
ically, the transfer functions O(s)/d(s) and 1y,(s)/o(s), which are the missile dynamics 
blocks in the autopilot block diagrams of Figure | on page | and Figure 2 on page 2, 
must be obtained. 

In order to determine these transfer functions, wind tunnel data and three sets of 
flight condition data are provided later m this chapter. The procedures for processing 
this data are also presented later mn this chapter. 


Prior to developing the nussile parameters, an explanation of the notation, termi- 


nology, and coordinate system used ts 1 order. [Ref ts py Ma2iiseler to Pictraues 





Figure 3.  Miissile Notation 


The notation in Figure 3 represents the following: 
eG is the center of gravity of the missile. 
tes the aerodynamic center of pressure of the mussile. 
ees the velocity of the missile. 
SRL is the spatial reference line. a reference fixed in space. 
F, is the lift force. It is perpendicular to the velocity vector. 
W is the weight of the missile. 
Mf, is the pitching moment. It is taken as positive in the counter-clockwise direction. 


o 1s the angle of attack. It is the angle taken from the velocity vector to the longi- 
tudinal missile axis. 


7 is the angle between the SRL and the longitudinal missile axis. 


6 is the fin deflection angle. It is taken from the longitudinal missile axis to the control 
fin axis. 


fp is the angle between the gravity vector and the SRL. 


Other notation and terminology will be discussed as it is introduced. It is important to 
remember that this problem considers only motion in the pitch plane but that similar 


methods may also apply to yaw motion. 


B. FLIGHT CONDITION DATA 

The three sets of flight condition data given in Table | on page 6 establish the de- 
sign points around which the autopiloted missile will be examined. The flight condition 
data will be used, in conjunction with wind tunnel data presented in the next section. to 
establish the missile parameters and autopilot gains necessary to develop and conduct 
simulations of the autopiloted missile. 

Some manipulation of the flight condition data is required in order for it to be used. 


These conversions and calculations will be made as necessary and explained when used. 


Table 1. FLIGHT CONDITIONS AT THREE DE- 
SIGN POINTS 


DESIGN Zoi wie 


Time (referenced to launch) = 9.5 s 
Alt = Altitude = 4.993 km 
vi Velocity = 1112.4 ms 
M Mass = 450 kg 


= Dynamic Pressure = 
I = Moment of Inertia = 
Pate Se (Orn 

CP 1s) 265 emitter mesc 


455642 N/m? 
TO Ken 


= 1416829 SS bee 
I = 693 kg om 
POWCE OND 
CP 1s 265 em ait o1moce 


eed 
ot 


] 
Hey ud 


AD tt ole 


Power Or 
CP 1s 266 cm aft of nose 





C. WIND TUNNEL DATA 

The wind tunnel data presented in this section will be instrumental in determining 
mussile parameters. The data. shown in Table 10 on page 37 through Table 13 on page 
40, are extracted from the complete wind tunnel data set for the missile used in this re- 
search paper. These data. as well as the three sets of flight condition data contained in 
Table 1, were provided by the U.S. Army’s Missile and Space Intelligence Center 
(STG): 

The first two sets of wind tunnel data are the axial force coefficients. These are 


given in Table 10 and Table 11 for the Power On and Power Off cases. 








The third set of wind tunnel data is the normal force coefficients. These data in 
Table 12 on page 39 will be used, along with the axial force coefficient data to calculate 
the lift and drag of the missile at the design points. 

The final set of wind tunnel data is the pitching moment coefficients. These data. 


in Table 13 on page 40 will be used to determine the pitching dynamics of the missile. 


D. WIND TUNNEL DATA INTERPOLATION 

To enter the wind tunnel data at the desired design point, the Mach number must 
be calculated from the flight condition data. Using given missile velocity and altitude, 
with temperature obtained from gas tables and graphs found in [Refs. 2,3], the Mach 
number for each design point 1s calculated by gas law equation techniques. These cal- 
culations are not included here. The Mach number at each design point is given in 


Table 2. 


Table 2. MACH NUMBERS AT DESIGN POINTS 


a a 






tJ 


The availability of the Mach number at each design point allows entry of the wind 
tunnel data and interpolation of the various coefficients for each design point. The in- 
terpolation of the wind tunnel data at each design point is accomplished using a simple 
linear interpolation scheme. The interpolation provides the coefficients at each design 
point. These interpolated coefficients are given in Table 14 on page 41 through 


Table 17 on page 43. 


E. CALCULATION OF LIFT AND DRAG COEFFICIENTS 
Once the axial and normal force coefficients at the design points have been deter- 
mined the lift and drag coefficients can be calculated. The lift force, F,, and drag force, 


F, , are related to the normal force, Fy, and the axial force, F,, by the trigonometric 


relationships 
ry, =f ycosa — sine (2p 
and 


Vy = Fxsina + F cosa. (oe 


These relatronships are depicted in Figure 4 below. 





Figure 4. Force Relationships 
Since the coefficients are dimensionless representations of the forces, the same re- 
lauionships apply. Specifically, 
C, = Cacosa — C, Sino (23) 
and 


Cy = Casing + C,cosa (2.4) 


Perc is tne miitcoelticient,@,us the drag coeflicient, C, is the axial force cocfficient, 
and C, is the normal force coefficient. [Ref. 1: p. 1344] 

The application of these equations to the interpolated wind tunnel data in Appendix 
B renders lift and drag coefficients at the design points. These lift and drag coefficients 


are tabulated in Appendix C. 


F. ANALYSIS OF WIND TUNNEL DATA 

There are four missile parameters which are of interest in the development of the 
mussile dvnamics transfer functions. These parameters are: F,,, which is the lift force per 
mimelewot attack, £,,, which is the lift force per fin deflection angle, :\7,,, which is the 
moment about the center of gravity per angle of attack, and 4/,,, which is the moment 
about the center of gravity per fin deflection angle. 

These parameters can be determined from their corresponding coefficients and a few 
physical dimensions of the missile. The corresponding coefficients are discussed in the 
paragraphs which follow. 

C,, is the lift coefficient per angle of attack. It is determined bv graphing the lift 
coefficients against the angle of attack and approximating the slope of the nearlv linear 
graph. 

C,,1s the lift coefficient per fin deflection angle. It is determined by graphing the lift 
coefficients against the fin deflection angle and approximating the slope. 

Similarly, Cy, and C,,,; are the pitching moment coefficients per angle of attack and 
fin deflection angle. respectively. They too are obtained by graphing the coefficients 
against the angles and approximating the slope of the resulting graph. 

The graphs of all the coefficients versus corresponding angles are contained in Ap- 
Bera. . The graplis are made for each set of coefficients at each of the three design 


mois | ne cociiicients C,,, C,,, Cy, , and C,,, are tabulated in Table 3 on page 10. 


Table 3. COEFFICIENTS DUE TO ANGLE OF AlLIAGRAND 
FIN DEFLECTION ANGLE 


DESIGN DESIGN DESIGN 
PO aes ie 





G. CALCULATION OF MISSILE PARAMETERS 

As mentioned in the preceding section. the parameters of interest are /,,, P77 
and M,;. With knowledge of the lift and moment coefficients previously obtained, and 
knowledge of some basic physical dimensions of the mussile, these parameters are deter- 
nuned. The necessary iissile dimensions are shown in Figure 5. These missile dimen- 


sions were provided by the Missile and Space Intelligence Center. 


SIDE VIEW 


5.3 m 


Se 0.38 m 
END VIEW 


Figure 5. Missile Dimensions 
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The equations used to determine the missile parameters from the coefficients are: 


OE Seis (2.5) 
Fig = Crs, 226) 
MW Onl ee (2.7) 
and 
Mpg = CyysqAc (2.8) 


where g is the dynamic pressure, provided in the flight conditions, A is the missile char- 


acteristic area, and c 1s the characteristic length. The characteristic area, A, 1s taken as 





(2.9) 


perere cd 1S the missile diameter. 38 cm. as given in Figure 19 above. The characteristic 
length 1s taken as the missile diameter. [Ref. 1: p. 134] 

Applving the above equations to the Known coefficients produces the missile pa- 
rameters at the three design points. Use of these parameters in later chapters requires 
that thev be given in force or moment per radian rather than per degree. This conver- 
sion 1s made mulitplving each parameter bv [8O/z . These converted parameters are 
tabulated in Table 4 on page 12. 

Knowledge of these missile parameters will enable the derivation of the missile dv- 


namics transfer functions in the following chapter. 


DI 


Table 4. MISSILE PARAMETERS AT DESIGN POINTS 
a 
ee FRG SS 238330 99695 N 
ee 


















Hl. DETERMINATION OF TRANSFER FUNCTIONS 


A. OVERVIEW 

In this chapter, the missile parameters obtained in the previous chapter are applied 
to the equations of motion which govern the flight of the missile. From this, the missile 
dvnamics transfer functions, 0(s)/6(s) and 77,,(s)/d(s), are derived. By then inserting these 
two missile dynamics transfer functions into the two autopilot transfer functions, the 
complete autopilot transfer functions can be expressed in terms of missile parameters 


and feedback gains. 


B. MISSILE DYNAMICS TRANSFER FUNCTIONS 


Consider the forces and moments acting on a missile while in flight. These are de- 


x 


picted in Figure 3 on page 4. Summing the forces in the transverse direction, that 1s 


Peroenaicular to the velocity vector, vields 


») Fre Wyn = LF, — WW sin(f + 6-2). (oa) 
Since 
Fi, =Ff,,4+F, (0 + a), (S22) 
the substitution is made yielding 
Ul se (Or Oh ge NO ar ep 1 belles ae ee E (Sea) 
Dividing through by the mass, m, yields the lateral acceleration, 7, . 


le ie WV 
Ny =r a + (5 + a) — — sin(B + 6 ~ a) (3.4) 


The weight of the missile is neglected at this point for simplicity. The missile’s 
weight would affect the lateral acceleration of the missile most when the angle 
(8B +@—«a) is 90 or 270 degrees. At this angle the contribution to lateral acceleration 


would be plus or minus 1 G. This effect would be offset by fin deflection and angle of 


attack trim angles of much less than one degree at any of the three design points. 


Neglecting missile weight leaves the approximation 





ge Sane 5 (Sear 2 Aaa - 
|= sa © es ee ee n> (335) 


which is used to determine 7,, at each design point in terms of ~ and 0. 
In a similar manner, summing moments about the center of gravity of the missile 


with the counter-clockwise direction being referenced as the positive direction yields 


> Meg = 16 — Mp (3.6) 


where J is the missile’s moment of inertia, @ is the missile’s angular acceleration in the 
pitch plane, and A/, is the pitching moment. The pitching moment, \/,, 1s given by the 


equation 
Mp = 6M, + OMy, + oMp, + (a + 5) Mpg. (3.7) 
Since « = 6 and & = @ the substitution is made leaving 
Mp = (Mp, + M,,) + o( Mp, + Mp5) + 5(Mps). (3.8) 
Combining and rearranging the above equations vields 


Mt Mp) (Mog + Mos) Mp 
“ —-———— oi FE 


: ; Po. (3.9) 


In Equation (3.9), the term 


(Mp; + Mya) 
[ 


is the aerodynamic damping term. It is bevond the scope of this research paper to 
properly obtain this term. It is, therefore, assumed that the missile is designed with a 
damping coefficient of €=0.5. The aerodynamic damping term will be different nu- 


merically at each of the three design points and will be referred to hereafter as D. It will 


be calculated later in this chapter. (Ref. 4] 


Substituting D into Equation (3.9) vields 


Leone. ane 
er oe (3.10) 


Using LaPlace transforms and assuming Zero initial conditions this equation transforms 


into 


(Mp, te M ps) = Mops es 














s*a(s) — Dsa(s) - eres cae a(S) — 7 6(s). (3a) 
Rearranging, this becomes 
Mops 
eS : | (3.12) 
O(S) 5 lls ae Uy 
Spe 
I 
Multiplying both sides by s renders 
Mops 
se.(s) 7 “(s) aes 
i = a a Le = nes) 
6(s) ; (Mp, + Mps) —— 6(s) 
Spe 
if 
Since & = 8, 
y 0 
2S) = oi) ———————— , (3.14) 
d(s) — A(s) 5 (Mp, + Mos) 
— Ds—- 7 


This transfer function 1s depicted in the signal flow graph of Figure 6 on page 16. 


—" 
Ler 





Figure 6. Pitch Transfer Function Signal Flow Graph 


The aerodynamic damping term, D, is calculated for each of the three design points. 
The denominator of the 6(s)/d(s) transfer function is the characteristic equation of that 


system. It 1s in the quadratic form of s?+ s2¢w, + w?. Equating coefficients yields 
—-D=2o, (3.15) 


and 


NT pa + Af 
;. | (3.16) 
I 
Using the missile parameters obtained in the previous chapter and the assumed value 
of 0.5 for ¢, the damping term, D, is calculated. The values for D are given in Table 5 


on page 17. [Ref. 5: pp. 106-107] 
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Table 3. AERODYNAMIC DAMPING COEFFI- 
CUEINS 


DES GN OLS 1 Denim COMP Cite | 





The other nussile dynamics transfer function needed is the lateral acceleration 
transfer function, 7,,(s)/d(s). Combining the signal flow graph of Figure 6 with 


Equation (3.5) results in the signal flow graph of Figure 7. 





Figure 7. Acceleration Transfer Function Signal Flow Graph 


The transfer function associated with this signal flow graph 1s given as 


Mps (iat Firs) 
te ttt (3.17) 
, Chelan) a | 





The transfer functions for the missile dynamics are thus obtained. They will be used 


in the following section to determine the autopilot transfer functions. 


C. AUTOPILOT TRANSFER FUNCTIONS 

By inserting the two missile dynamics transfer functions into each of the two 
autopilot configuration block diagrams presented in Chapter One, and performing a few 
basic block diagram manipulations, the closed loop autopilot transfer functions are ob- 
tained. In both cases. the desired autopilot transfer function is 74,(s)/y-(s). For Autopilot 


Configuration One the transfer function is 








N af(S) 
1 (5) & o(s) (3.18) 
i1c(s) 5 ; A(s) . Q(s) : | | 
S + s(Kpt+ Rey) ae oie See a Ky 


Expanding the missile dynamics transfer functions in the autopilot transfer function and 


performing extensive algebraic manipulations vields 








nus) Ny 1S? ae AS tela 
G,{s) = Pans. = (3.19) 
TOD cial 6 Maa gs ara ah, 
where 
Fig 
imma ei f 
Fis 
ny == D m °? 
— Mpg (Frnt Fis) Fig (Mpa + Mpg) 
i ] 7 asi ] 
d,=Ap—D, 
, Mos smlans P (Mp, + Mpg) 
d\, = Kgp 7 aE iil an pD— i ; 
a) am AsrKpp ] ae Kr ] a K4D mt? 
and 


1S 


r ar Mn ] 


The transfer function for Autopilot Configuration Two is obtained in a similar 


manner. For the second configuration 














(Ss 
7 (Kos + Ky) atl ) 
144(S) 7 O(s) (3.20) 
TO) (Ket Kp oo) + KylKys + Ky) 
/ aa BR 6(s) Ax**2 } 6(s) 
This expands into 
H (5S) figgS” hates: + ios +n 
G,(s) = M == (3.21) 
NCIS) §° + dy + dos” + dyos + dy 
menere 
p 
aD) Ky + 
> Fis , Fis 
io = ar ari 
- pipe elie) Vigieenlelt, belies) 2 meen, 
i ann ee eens ce eee ) Da, 
pee ee Oe (UE a) 
N32 = Ki mr ort i 
' 6 aoe 
d,=Krp—D+K,fA, m ’ 
- Mp; Mp, ae Mp; ye ‘ Fy s Aer Fis 
diz = Karp > > OAD iD t+ Kai aa 
Ree raar, 5) Fi, (ip, + Mz) Fis (Mp, + Mops) 
dy = KK, a mM > iy ] )- K,K,D ie ] ’ 
and 
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These closed loop autopilot transfer functions will be used in the next section for 


system pole placement and solution of the feedback and feed forward gains. 


D. POLE PLACEMENT AND GAIN DETERMINATION 

As stated in Chapter One of this research paper, the desired autopilot response for 
both configurations is w, = 10 with a damping coefficient €=0.5. These design spec- 
ifications dictate a complex pair of poles in the characteristic equation whose real part 
is -5.0. Since the autopilots’ characteristic equations are fourth order, two more poles 


are required to put them in the form 
GES Wi), oe ese aoa). (323; 


In order for the complex pair placed by design specifications to be dominant. the other 
two poles must le to the left of them in the s-plane. The two additional poles are se- 
lected for simplicity as real and equal and placed at -20.0. The characteristic equation 


for both autopilot configurations 1s 
CE =(s + 20)(s + 20)(s* + 10s + 100) 


= 5° + 505° + 900s” + 8000s + 40000. (3.23) 


The denominators of the closed loop transfer functions, which are given by 
Equations (3.19) and (3.21). are the autopilot characteristic equations. The transfer 
function gains are solved for by first equating the coefficients of characteristic equation 
of above to those of Equations (3.19) and (3.21), then solving simultaneous equations 
for the gains. The calculated gains for Autopilot Configuration One and Autopilot 
Configuration Two are tabulated in Table 6 and Table 7, respectively. The determi- 
nation of the autopilot gains makes it possible to simulate the response of the missile. 


These simulations will be conducted in the next chapter. [Refs. 6,7] 


Table 6. AUTOPILOT GAINS FOR CONFIGURATION ONE 
EES DES DESIGN PLSIGN 
a 





Table 7. AUTOPILOT GAINS FOR CONFIGURATION TWO 


DESIGN Wiss iG DESIGN 
POINT 1 POINT 2 POEUN 3 





It 1s interesting to note that the values of K, and A, for Configuration Two, depend 
on the selected value of A,. For purposes of simulation in the next chapter, A, 1s chosen 
as unity, This results in K,A, being equal to A, and 4,4, being equal to A, for the sim- 


mlation. 


IV. AUTOPILOT SIMULATION AND ANALYSIS 


A. AUTOPILOT SIMULATION 

In order to compare the responses of the two autopilot configurations, a simulation 
program was developed using IBM's Dynanuic Simulation Language. [Ref. §] 

. The program simulates the step response of the autopilot transfer functions for both 
configurations, at all three design points. The simulation was developed from Equations 
(3.19) and (3.21), which describe the two autopilot transfer functions, and uses the 
feedback and feed forward gains tabulated in Table 6 and Table 7 on page 21. The 
program also incorporates a saturating limiter which limits the deflection of the control 
fins to + 20°. A copv of the program listing 1s at Appendix E. 

The program was run on the IBM 3360 mainframe computer system at the Naval 
Postgraduate School utilizing a Tektronix 618 Monitor and Tektronix 4631 Hard Copy 
Lnit for graphic output. 

This chapter contains the graphical results of these simulations. Figure 8 on page 
23 through Figure 10 on page 25 depict the step response of Autopilot Configuration 
One at each of the three design points. Figure 11 on page 26 shows the parametae 
display of the step response of Autopilot Configuration One at all three design points. 
Figure 12 on page 27 through Figure 14 on page 29 depict the step response of 
Autopilot Configuration Two at each of the three design points. Figure 15 on page 30 
displays parametrically the step response of Autopilot Configuration Two at all three 
design points. The parametric comparisons of the two autopilot configurations at each 
design point are shown in Figure !6 on page 31 through Figure 18 on page 33. The 


neXt section contains discussion and analvsis of the graphical output. 
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Figure 8. 
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Step Response of Autopilot Configuration One at Design Point One 
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Figure 9. 
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Step Response of Autopilot Configuration One at Design Point Two 


Figure 10. 
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Step Response of Autopilot Configuration One at Design Point Three 
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Figure 11. Parametric Comparison of the Step Response of Autopilot Configuration 


One at Three Design Points 


1.8 2.9 


1.6 


1.2 


0.8 


0.6 


0.2 


0.9 


213020 


Figure 12. Step Response of Autopilot Configuration Two at Design Point One 
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Step Response of Autopilot Configuration Two at Design Point Two 
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Figure 14. Step Response of Autopilot Configuration Two at Design Point Three 
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Figure 15. Parametric Comparison of the Step Response of Autopilot Configuration 


Two at Three Design Points 
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Figure 16. 


Parametric Comparison of the Step Response of Autopilot Configura- 
tions One and Two at Design Point One 
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Figure 17. Parametric Comparison of the Step Response of Autopilot Configura- 


tions One and Two at Design Point Two 
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Figure 18. Parametric Comparison of the Step Response of Autopilot Configura- 


tions One and Two at Design Point Three 
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B. ANALYSIS OF SIMULATION RESULTS 

In the analvsis of the sunulation results, there are two characteristics of the step re- 
sponse which are of interest. Thev are the rise time, z, , and the setthing time. ¢,. 

The rise time is simply defined as the tme required for the systems output response 
to go from 10 percent to 90 percent of its final value. This can be measured directlv off 
of the graphical simulation results. [Ref. 9: p. 40] 


- The settling time, ¢, , is defined as four time constants, or 





(4.1) 


The value of the maximum overshoot, .J/,,, and the time at which the maximum over- 


SlOOtOcCUIS. ¢ 


, » can be measured directly off the graphical output. The maximum 


overshoot 1s defined as 


Y 


Rs 


M,, = 1+ exp( == (4.2) 


and the time at which it occurs 


=< : (4.3) 


On l—t 

Measurements of \f,, and 1, are made from the graphical output of Figure § on page 
23 through Figure 10 on page 25 for Configuration One and from Figure !2 on page 
2/through Figure 14 on page 29 for Configuration Two. Once JZ, and 1, are measured, 
¢ and w, are obtained by solving simultaneously Equations (4.2) and (4.3). The values 
for ¢ and w, are then used in Equation (4.1) to obtain the settling time, z,. The response 


characteristics JZ, 1, 6, w, , #, and /, are tabulated in Table 8 and Table 9 on page 35. 
[Ref. 5: pp. 108-109] 
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Table 8. RESPONSE CHARACTERISTICS FOR AUTOPILOT 
CONFIGURATION ONE 


DE clGes DEeslGy, DESIGN 
POR we FOUN: POENT 3 















Table 9. RESPONSE CHARACTERISTICS FOR AUTOPILOT 
CONFIGURATION TWO 


IDES CON De SGN DEstiGe 
HOON TE OAS AT EOIN 
ae 
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V. CONCLUSIONS 


The step responses of the two autopilot configurations were compared at three dif- 
ferent design points. Both configurations met the desired response specifications of 
¢~0.5 and w,~10. The actual values of € and w, are given in Table 8 and Table 9 on 
page 35 for Configurations One and Two respectively. 

At Design Point One, Configuration One’s observed rise time is nine percent less 
than that of Configuration Two, and Configuration One’s observed settling time is 
twenty-three percent less than that of Configuration Two. The two configurations are 
compared graphically in Figure 16 on page 31. 

At Design Point Two. Configuration One’s observed rise time is equal to that of 
Configuration Two. and Configuration One’s observed settling time is nineteen percent 
less than that of Configuration Two. The graphical comparison is shown in Figure 17 
Oh pace... 

At Design Point Three, Configuration One’s observed rise time 1s sixtv-four percent 
more than that of Configuration Two, and Configuration One’s observed settling time 
is thirty percent less than that of Configuration Two. This comparison is shown graph- 
ically in Figure 18 on page 353. 

The apparent trend is that Configuration One’s rise time increases slightlv as the 
flight of the missile continues but that its settling time decreases slightly. This trend is 
graphically depicted in Figure 11 on page 26. Configuration Two’s rise time decreases 
sigmificantly as the flight of the missile continues but its settling time remains nearly 
constant. This trend is shown graphically in Figure 15 on page 30. 

The shorter rise time, later in flight, of Configuration Two. is viewed as a significant 
performance advantage over Configuration One. It is clear that a shorter rise time 
means a faster response. This faster response is desirable late in flight because it is 
during this later part of the missile’s flight that the missile would be in its final phase of 
closing on and tracking a maneuvering target. For this reason, faster response, later in 


flight, is a significant advantage. 
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APPENDIX A. ORIGINAL WIND TUNNEL DATA 


Table 10. AXIAL FORCE COEFFICIENTS (POWER ON) 
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Table 11. AXIAL FORCE COEFFICIENTS (POWER OFF) 
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Table 12. NORMAL FORCE COEFFICIENTS 
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Table 13. PITCHING MOMENT COEFFICIENTS 


DELTA | Err | MERGE 2 CO |i @inecnagn acesvenigersn 
(degrees) 
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APPENDIX B.) INTERPOLATED WIND TUNNEL DATA 


Table 14. ANTAL FORCE COEFFICIENTS (POWER ON) 
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Table 15. ANTAL FORCE COEFFICIENTS (POWER OFF) 
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Table 17. PITCHING MOMENT COEFFICIENTS 
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APPENDIX C. LIFT AND DRAG COEFFICIENT DATA 


Table 18. LIFT COEFFICIENTS (POWER ON) 
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Table 19. LIFT COEFFICIENTS (POWER OFF) 


DELTA ALPHA DESKS. DESIG. DESIG. 
(degrees) (degrees) POSS! LOIN T 2 ROUNT > 
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Table 20. DRAG COEFFICIENTS (POWER ON) 
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Table 21. DRAG COEFFICIENTS GMa ER OFF) 
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APPENDIX D. WIND TUNNEL GRAPHS 
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Figure 19. Lift Coefficient versus Angle of Attack at Design Point One 
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Figure 20. Lift Coefficient versus Fin Deflection Angle at Design Point One 
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Figure 21. 
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Pitching Moment Coefficient versus Angle of Attack at Design Point 
One 
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Figure 22. Pitching Moment Coefficient versus Fin Deflection Angle at Design 
Point One 
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Lift Coefficient versus Angle of Attack at Design Point Two 
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Figure 24. 
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Lift Coefficient versus Fin Deflection Angle at Design Point Two 
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Figure 25. 
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Figure 26. Pitching Moment Coefficient versus Fin Deflection Angle at Design 
Point Two 
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Figure 27. 
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Lift Coefficient versus Angle of Attack at Design Point Three 
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Figure 28. Lift Coefficient versus Fin Deflection Angle at Design Point Three 
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Pitching Moment Coefficient versus Angle of Attack at Design Point 
Three 
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Figure 30. Pitching Moment Coefficient versus Fin Deflection Angle at Design 
Point Three 
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APPENDIX E. SIMULATION COMPUTER PROGRAM 


: a we'selseaclactan ses seve «sls rs wlan sols ulaeton s's s's x's =": wan selaselacts 
ves ves 4h Martel tarde ote its seveses'esies 


setede THESIS RESEARCH Ae 
seteds DSL SIMULATIONS FOR TWO Jeter 


Bor AUTOPILOT CONFIGURATIONS verse 
oe AT THREE DESIGN POINTS weaere 


wereve KENNETH E. COCKERHAM Fereve 
Waev CAPTAIN, UNITED STATES ARMY Cee 


meas NAVAL POSTGRADUATE SCHOOL wrere 


7. risks Wvavevevevevevaev ee Tee TONS re ie ier ivwielwiltwiwiwitiwiwiwi wis 


BV CN ON ON BN ON ON BREN ON SN ON ON ON ON ON Se ae Oe ON ED Ev ON OD NOD rig rier erie iwi eri ie iy A OV OV OV ON ON ON AN ON AN ON AN 


FIXED INPUT 
oN 1 
Youle stes'e Jee sles touJe Seorlerten te ves “DES T GN PO INT ONE* ale steslesfectociest afoslecteslonle Jerse Sestestevlovle Je 


vetevededesededs se veved sees MT SST LE PARAMETERS wvstede cost vostvevedevoredodedeceveve 
PARAM I = 720. O, M = 450.0, DAMP = -22.4, FLA = 799276.0, FLD = 88808.0 
ee MPA = - 168736. 0, MPD r= 92191365720 


sseveseseirsest FEEDBACK GAINS FOR CONFIGURATION ONE? stsescererssese 
ae KAl = -0.094, KF1 = 27.6, KBR1 = 0.752, KST = 26.96 


were FEEDBACK GAINS FOR CONFIGURATION TWO *eressteessorse 
PARAM KAZ = 1.0, KF2 = 20°71, KBRZ =e -O iS 


ae “*FEED FORWARD GAINS FOR CONFIGURATION TWO*ssvedvecers 
PARAM Se = es pee ep = ine oe 


= v8 ae aes 5 TONS aaa ae GS aN EN OV ON GV EVEN EV EV EVEN DREN ED OVEN CRON ON ON OD 


velevederededeteteee AUTOPILOT CONFIGURATION ONEvsitteveveretetevevererere: 


KAl * STEP(0. 0) 
Rivcakige= sa 
PNIGEL (O10. an) 
X2 - X11 - X13 
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X4 INTGRECO. 0 ,X3) 


K5 = X4 

DELTAl = LIMIT(-0. 35,0. 35,X5) 
X6 _ = DELTA1 * MPD/TI 
ALFDD1 = X6 + X7 + X8 

ALFD1 = INTGRL(0.0,ALFDD1) 
ALF 1 = INTGRL(0. 0,ALFD1) 
X7 = ALFD1 * DAMP 

X8 = ALF1 * (MPA+MPD)/I 
X9 = DELTA1 * FLD/M 

X10 = ALF1 * (FLA+FLD)/M 
ACCEL! = X9 + X10 

mii = KBR1 * ALFD1 

X12 = KST * X11 

X13 = Om ee 

X14 = KAl * ACCEL1 


ve 
waedekdededevedevee AUTOPILOT CONFIGURATION TWO wrrvedswvssdededededed sve 


Xe 


R2 = KA2 * STEP(0. 0) 

Y1 = R2 - Y16 

Y2 = Kil * Yi 

Y3 = INTGRL(0. 0, Y2) 

Y4 SRD ee 

Y5 = Y4 + ¥3 

Y6 = Y5 - Y14 - Y15 

Y7 = INTGRL(0.0,Y6) 

Y8 = Y7 

DELTA2 = LIMIT(-0. 35,0. 35,Y8) 
y9 = DELTA2 * MPD/I 
ALFDD2 = Y9 + Y10 + Y11 
ALFD2 = INTGRL(0.0,ALFDD2) 
ALF2 = INTGRL(0. 0,ALFD2) 
Y10 = ALF2 * (MPA+MPD)/I 
vata = ALFD2 * DAMP 

fale = ALF2 * (FLA+FLD)/M 
wile = DELTA2 * FLD/M 
ACCEL2 = ¥12 + Y13 

Y14 = ALFD2 * KBR2 

Y15 =) Yen 

Y16 = KA2 * ACCEL2 


oc 
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CONTROL FINTIM =2.0, DELT = 0.0001 
ve 

SAVE 0.001, ACCEL1,ACCEL2 

END 


se 


we 
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bh 
aaah me 

vedevederevedevee Ie DESTGN POINT TWO %sttedtsededeckvedsvededovede ves pee vw 
Tevedededededs tevededtekeeM T SST LE PARAMETERS **%%* Sovevedevedesevesewetss 
PARAM I = 693. M S2370 -2Dshr = eS. FLA = Daa teaion FLD = 19080. 0 
PARAM MPA = -43545.0, MPD = -54431.0 
wisdevedereteeFEEDBACK GAINS FOR CONFIGURATION ONE exer 
PARAM KAl = -0.845, KF1 = 38.1, KBR1 = -4.44, KST = 8.98 
waswweweues FEEDBACK GAINS FOR CONFIGURATION TWOrrieveredesesesest 
PARAM KA2 = 1.0, KF2 = 41.05, KBR2 = -7. 62 
wersewese FEED FORWARD GAINS FOR CONFIGURATION TW Over veces aece 
PARAM kl = -0.845, K2 = -0.057 
Bot orate ie oe oe ae) ei) eae econ Sesetess Teme sere reese reset aC Ue CU ee TCO Teen nae 
END 
a 
chy 
ry 
pee 3 
sees Lar ierlscharlae hockii eee POINT THREE “vere sttededevededededestdevede ss ves 
sevedevevevedededs MISSILE PARAMETERS * eeu Pie Liuirbiri arias oe Sirs 
PARAM I = "687, M = 365, DAMP = -9.1, FLA = = . 99695. 0, FLD = 10657.0 
eae os ~ 26986. OO} MER =—-2767 0 


smear seacee “FEEDBACK GAINS FOR CONFIGURATION ONE? ?8ecsosesrecese 
PARAM Kal = =3706, KE i= 40 eel eee “XST = 10.18 
Wasser wae FEEDBACK GAINS FOR CONFIGURATION TWOsnirnieeeiere 


PARAM KA2 = 1:0.° KE2 = S11 lye eRe = 

wesw REED FORWARD GAINS FOR CONFIGURATION TWO**%eseesr8 
PARAM Kl = ~3. ee K2 = oe Soc 
Sevevevedevevevevevesesevedevesevedeveslevosevesesese ce seve Fete Tes e Te Te TOT CTC TEE TO FETE UE OT CTC IETS 


END 
GRAPH (Gl ,DE=TEK618) TIME, ACCEL1(RU=1,LO=-0. 2,SC=0. 2) 
LABEL (G1) STEP RESPONSE, CONFIGURATION 1, DESIGN POINT 1 


GRAPH (G2 ,DE=TEK618) TIME, ACCEL1(RU=2) 
LABEL (G2) STEP RESPONSE, CONFIGURATION 1, DESIGN POINT 2 


GRAPH (G3 ,DE=TEK618) TIME, ACCEL1(RU=3) 
LABEL (G3) STEP RESPONSE, CONFIGURATION 1, DESIGN POINT 3 


GRAPH (G4 ,DE=TEK618) TIME, ACCEL1(RU=1, 2,3) 
LABEL (G4) STEP RESPONSE, CONFIGURATION 1, DESIGN) OlNtseiy 2) cs 


GRAPH (G5 ,DE=TEK618) TIME, ACCEL2(RU=1) 
Me (G5) STEP RESPONSE, CONFIGURATION 2, DESIGN POINT 1 


GRAPH (G6 ,DE=TEK618) TIME, ACCEL2( RU=2) 
LABEL (G6) STEP RESPONSE, CONFIGURATION 2, DESIGN POINT 2 
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LABEL (G7) STEP RESPONSE, CONFIGURATION 2, DESIGN POINT 3 
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LABEL (G8) STEP RESPONSE, CONFIGURATION 2,DESIGN POINTS 1, 2, & 3 
GRAPH (G9 ,DE=TEK618) TIME, ACCELI(RU=1 ,LO=-0. 2 ,SC=0. 2) ,ACCEL2( RU=1) 
LABEL (G9) STEP RESPONSE, CONFIGURATIONS 1 & 2,DESIGN POINT 1 

a: 

OeerH (G10 ,DE=TEK618) TIME, ACCEL1I( RU=2) , ACCEL2(RU=2 ) 

fee (G10) STEP RESPONSE, CONFIGURATIONS 1 & 2, DESIGN POINT 2 

iy 

GRAPH (G11 ,DE=TEK618) TIME, ACCEL1(RU=3) ,ACCEL2(RU=3) 

LABEL (G11) STEP RESPONSE, CONFIGURATIONS 1 & 2, DESIGN POINT 3 


Yn 
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END 
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